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ABSTRACT 


The rapid growth on technological frontiers and 
increasing competition in the market has adversely influenced the 
manufacturing plants turning out production in mass or batches of 
large volume. More and more companies are submitting to group 
technology and flexible manufacturing systems. The benifits 
accrued from a change over to group layout has made cellular 
manufacturing attractive to both academicians and practitioners. 
This study addresses the problem of cell formation in cellular 
manuf ac t ur i ng . 

Numerous heuristic and analytical methods have been 
developed and applied in literature, but very little effort has 
been put into comparing and selecting the best or appropriate 
method. Solution of practical problems without adequate 
generalizations and theoretical formulations is limitedly useful. 
When researchers and practitioners claim success with a certain 
method it is necessary to compare and evaluate result on absolute 
quantitative scale. 

The present work compares some selected algorithms on 
the basis of grouping efficiency, grouping efficiacy, measure of 
effectiveness, time, number of bottleneck machines and number of 


exceptional parts. 
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CHAPTER I 


INTRODUCTION 


1.1 Introduction 

Over the last two decades, group technology CGT) has emerged 
as an important scientific principle CKusiak [41]) in improving 
the productivity of manufacturing systems. Cellular manufacturing 
(CM) / GT can be called as one of the most important innovations 
of the last few years. More and more manufacturing industries 
involved with small lots and a variety of products are becoming 
interested in group technology, which is particularly applicable 
in the area of batch type manufacturing. 

Group technology has also been recognized as an essential 
element of the foundation for the successful development and 
implementation of computer aided manufacturing (CAM) through the 
application of the part family concept. Recognizing the 
contribution that GT can make to the implementation of Computer 
Integrated Manufacture (CIM), Merchant says, " the appropriate 
initial step, to lay a sound foundation for the gradual evolution 
of a factory to full computer control, is to institute group 
technology cellular organization*'. The application of -group 
technology does not depend upon the degree of automation, and may 
be applied in a totally automated or even in a manual production 
system (Choobineh [19]). The basis for implementation of cellular 
manufacturing systems is the allocation of parts or components 
with certain 'similarities’ into groups or families, and the 
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formation of machining cells which are capable of manufacturing 
one or more part families. Manufacturing control can then be 
exercised within or between cells, rather than at the individual 
machines, using a combination of computer controlled machining 
and handling devices. 

Group technology has probably had a greater impact on 
increasing manufacturing productivity ^ than any other 
manufacturing concept. This can be attributed partly to 
contributions made by cellular manufacturing concepts to other 
manufacturing technologies such as roboticised and flexible 
manufacturing systems. 

Group technology is a manufacturing philosophy in which 
similar parts are identified and grouped together to take the 
advantages of their similarities in manufacturing and design. The 
idea behind GT is to decompose a manufacturing system into 
subsystems. In addition to the simplification of management 
control through the creation of similar sub systems called cell, 
GT also leads to many other benifits. A cell is essentially a 
group of machines and a family of related components. For 
example, for a plant producing 10,000 different parts it may be 
able to group the vast majority of these parts into 50 or 60 
distinct families having similar design and/or manufacturing 
characteristics. This leads to similar processing of members and 
high manufacturing efficiencies. The cell is self contained and 
self managed, with performance control entrusted to the foreman 
who wields authority to change loading on machines, assigns 
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operators to machines or alter the sequence of operations that he 
deem fit for proper load balancing within cell so long as such 
decisions do not contradict the overall strategies laid down by 
top management. 

The part families are assigned to cells to satisfy some of 
the following objectives 

(1) Minimize total material handling cost which includes costs 
due to inter cellular moves and intra cellular moves. 

(2) Maximize within cell utilization of machines. 

(3) Minimize duplication of machines in different cells. 

(4) Minimize number of exceptional parts and inter cellular 
trips. 

With more effective design rationalization and manufacturing 
standardization, the efficiency associated with flow line 

production is sought while the flexibility of a job shop 

manufacturing system is maintained. These features are essential 
for a firm to remain competitive in the current manufacturing 
environment, in which more special orders are demanded, 
production life cycles are reduced, and competition is focussed 
on factors like delivery speed, quality, design flexibility, 
delivery reliability as well as price. Firms have to manufacture 
products in a larger product mix, smaller volume, increased part 

complexity and shorter production period. Group technology, with 

its promises to reduce lead time, material handling, set up time, 
expediting, work in process and finished goods inventories, was 
developed as an approach to address the difficulties faced by 
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today's manufacturing industry. More over the need for 
flexibility and Just in time processing capabilities is forcing 
many discrete parts manufacturers to consider reorganizing their 
facilities. 

The vast majority of modern engineering factories are 
required to manufacture many different products which are ordered 
irregularly and usually small quantities. The product range is 
usually wide and the company is forced to manufacture many 
thousands of different components. In such companies heaps of 
components are seen lying around the workshop floors. Components 
which are urgently required hold up the assembly stages of 
production. 

The inefficiency and disadvantages encountered for large 
variety and medium/small batch size production in the traditional 
manufacturing environment pose three important questions 

(1) What should be the type and number of facilities ? 

(2) What should be the layout of facilities ? 

(3) How should parts be scheduled and released for production ? 

The answers to the above three questions lie in the emerging 
concepts of group tephnology/cellular manufacturing system (CMS) 
and flexible manufacturing system. 

1.2 Implementation of group technology 

A gradual implementation of GT in stages in a factory is 
more commonly adopted when compared to in one step, because of 
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the inherent risks and disruptive effects to total 
re-organization. A GT programme could require 2-3 years to 
install. Generally, GT implementation team is created whose 
members are drawn from a variety of departments such as 
production, production planning and control and design. First a 
pilot line is established until a satisfactory long term solution 
is reached which suits the unique circumstances of the factory. 
At the same time existence of cell enables all members of the 
factory to become slowly accustomed to the new manner of work and 
re-organize the advantages. Additional cells are slowly 
introduced one by one till all components are manufactured within 
the cells. 

1.3 Scope of present work 

During the past decade, there has been tremendous interest 
in GT problem from practitioners as well as from academics. 
Numerous heuristic or analytical methods have been developed and 
applied, but very little effort has been put into selecting and 
comparing the best or appropriate method. 

Solution of practical problems without adequate 
generalizations and theoretical formulations is limitedly useful. 
When researchers and practitioners claim success with a certain 
method it is necessary to compare and evaluate result on absolute 
quantitative scale. 

Present work codes various algorithms and compares them on 
similar quantitative and qualitative scales. 
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1.4 Organization of the thesis 

Chapter I gives the introduction to the subject of group 
technology. The advantages and disadvantages of group technology 
are mentioned. The explanation of ’How to implement group 
technology’ is given. 

Chapter II gives the detailed literature survey. A brief 
observation on literature survey is made towards the end of the 
chapter. 

Chapter III explains various algorithms which are coded and 
compared in this work. 

Chapter IV describes various measures of comparisons on the 
basis of which various algorithms are compared. This chapter 
explains the concept of Grouping efficiancy, grouping efficiacy, 
measure of effectiveness, bottleneck machines and exceptional 
parts. 

Chapter V lists the observations and comparisons done while 
solving the illustrations with various methods. The explanation 
for the occurance of the observations is also given. Four 
illustrations, relevant tables and graphs are also given in this 
chapter. 

Chapter VI gives the conclusions and the suggestions for the 


further work. 
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Appendix A gives the calculations of the various methods 
explained in the chapter 3. It also shows the calculations for 
calculating grouping efficiancy and grouping efficiacy. 



CHAPTER II 


LITERATURE SURVEY 


In this chapter a brief description of some of the major 
efforts in the field of group technology is presented. Since 
group technology aims at dividing the production systems, most of 
the research is focused on dividing the machine and components 
into groups. To understand what we mean by dividing of production 
system, following example is given. 

Consider the machine part incidence matrix (Fig 2.1) in which 
= 1 indicates that component j visits machine i, and a.j is 
blank, otherwise. 
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The organization of description of major efforts in the 
field of group technology, is based on classification scheme. The 
efforts can be classified under approaches like 

(1) Rule of thumb 

(2) Classification and coding system 

(3) Flow analysis 

(4) Approaches using similarity coefficient 

(5) Cluster analysis 

(6) Mathematical programming approach 

(7) Heuristic approaches 

(8) Hybrid approaches 

(9) Fuzzy theoretic approaches 

(10) Expert System / Artificial Intelligence approach 

2.1 Rule of thumb 

Based on the production engineer’s local product knowledge 
of the parts and processes, some families of the components may 
be self evident. Families consist of components defined by name 
or functions e.g pump casings, shafts, and sleeves. 

2.2 Classification and coding system 

A classification and coding scheme sorts parts into 
different classes based on certain part characteristics (such as 
routings) and assigns a code (usually alphanumeric) such that 
parts having similar codes could be easily identified as parts 
having similar codes could be easily identified as part families. 
A number of commercial (for e.g BRISCH, CODE, MICLASS) as well as 
non-proprietary (for e.g OPITZ, 'KC - 1) coding schemes are 



currently available. Most schemes are useful for variety 
reduction, and design information retrieval, and rarely have the 
capability as stand alone procedures for finding part families 
for cellular manufacturing. This is due to the fact that they 
group parts which are similar in shape but differ in batch size 
or tolerances and, hence should be made on different types of 
machines. Additionally, they do not group parts which are similar 
in design features but are processed on the same set of machines 
(Burbidge [9]). There are a few exceptions. The MICLASS system 
developed by TNO in Holland considers both design as well as part 
routing features and can be used to form part families and 
machine cells for cellular manufacturing systems. 

2.3 Flow analysis 

The aim of flow analysis is that of ’finding the families of 
components and associated groups of machines for group layout, by 
a progressive analysis of the information contained in route 
card’ . 

2.3.1 Production flow analysis 

Production flow analysis (PFA) was first introduced by 
Burbidge [8]. It involves the systematic listing of the 
components in various ways in the expectation that groups of 
machines and components may be found by careful inspection. As 
Beer and de Witte [22] point out, the procedure requires a 
series of evaluations to be made by the designer, more or less 
calling upon his ability to recognize pattern. 



PFA consists of three levels of analysis. 

(a) Factory flow analysis 

It makes use of process route numbers, in order to get an 
overall picture of the present state of material flows. Machines 
are divided into departments, and each department is given a 
number. The process route number of a component is defined as the 
sequence of the numbers of the departments visited. A flow chart 
showing the interaction of various departments based on these 
numbers is then drawn. 

(b) Group analysis 

It attempts at effective formulation of groups of machines 
and their families of related components. The 'best" method 
called "Nuclear Synthesis’ is based on selecting machines used by 
a few components as starting points for various cells, or nuclei. 
The next machine is allocated on the basis that it has the 
smallest number of components left unassigned to a group once 
"Nuclear Synthesis* is completed, these nuclei are modified and 
combined, until the required number of groups is formed. 

(c) Line analysis 

It aims at creating a pure line-layout between machines 
inside the groups. 

2.3.2 Component flow analysis 

Component flow analysis (CFA) was first used by El-Essawy 
[26]. CFA employs three stages of analysis. 



FIRST STAGE : The objective is to consider the total component 
mix of the company and to identify and sort components into 
categories according to their manufacturing requirements. The 
components are sorted in order of machine requirements and sorted 
list are printed in two ways, firstly in the order of the number 
of machines required and secondly in the order of the smallest 
machine numbers involved. 

SECOND STAGE : The aim is to obtain groupings of the machines 
using the lists of sorted components and taking into account 
various local constraints. Rough groups are formed, to which 
other machines and components are successively added. 

THIRD STAGE : It involves a detailed analysis of the loadings 
and flow pattern of the cells with appropriate adjustments to 
ensure that an acceptable design is achieved. 

PFA differs from CFA in following ways. PFA first partitions 
the problem whereas CFA does not. The manner in which the cells 
were built up is also different in two methods. 

2.4 Approaches using Similarity coefficients 

The basis of this method is to measure the 'similarity’ 
between each pair of machines and then to group the machines into 
families based on their similarity measurements. In most cases, 
the similarity measurement used is the coefficient of Jaccard 
which is defined for any pair of machines as ' the number of 
components which visit both machines, divided by the number of 
components which visit at least one of the machines' . The 
components having similar coefficients are grouped together to 



form a cell. 


Of the several similarity coefficients selection of a 
particular similarity coefficient depends upon its application. 


M 

A 
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ASC. . 
ij 


SC 

"iJ 

TNC^ 
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: Total no. of machines in the system. 

: Number of machines used by both components. 

: Number of machines used by one of the component. 

: Number of machines used by the other component. 

: Number of machines not used by either of the 
components. 

: Additive type of similarity coefficient between 
i th and j th machines. 

: Similarity Coefficient. 

: Number of common components using both machines i and j. 

: Total number of components using i th machine. 

: Similarity coefficient of product type based on 
total number of common components processed by 
each of the machines i and j. 

: Total number of common components processed by 
i th machine. 


Similarity coefficient can be defined as 

SC = A / (A + B + C) (Waghodekar and Sahu [68]) (2.1) 

ASC. . = e. . / (TNC,+ TNC - e ) (2.2) 

ij ij 1 J iJ 

(Rajagopalan and Batra [54] ) 


Other similarity coefficients can be defined as 



(a) sc = A / M 


(2.3) 


(b) SC = A / (B + C) (2.4) 

(c) SC = 2A / (2A + B + C) (2.5) 

(d) SC = B + C / M (2.6) 

(e) SC = (A + D) / (B + C) (2.7) 


Other similarity coefficients can be defined as 

(a) SC = (e^j)^/ (TNC^» TNC^) (2.8) 

( if Ej = 0, assume 

(b) . (e^.)2 ^ ^ 

McAuley [47] uses single linkage cluster analysis. This 
method first clusters together those machines mutually related 
with the highest possible similarity coefficient, then it 
successively lowers the level of admission by steps of 
predetermined equal magnitude. The admission of a machine or 
groups of machines into another group is by a criteria of single 
linkage. The main disadvantage of this method is that while two 
clusters may be linked by this technique on the basis of a single 
bond, many of the members of the two clusters may be quite far 
removed from each other in terms of similarity. 

Rajagopalan and Batra [54] developed .a graph theoretic 
method which uses cliques of the machine-graph as a means of 
classification. The vertices of this graph are the machines and 
the arcs are the Jaccard similarity coefficients and a clique is 
a maximal collection of vertices, every pair of which is 
connected by an edge of the graph. The main disadvantage of this 



approach is that because of the high density of the graph, a very 
large number of cliques is usually involved and many of the 
cliques are not vertex disjointed. To reduce the number of groups 
and to incorporate the machines which are not included in the 
cliques, graph partitioning is used and it is at this stage that 
the allocation of components, in accordance with a number of 
heuristic rules, is also carried out. 

2.5 Cluster analysis 

Cluster analysis is the widely used technique for cell 
formation in group technology. These methods re-arrange the rows 
and columns of the machine component incidence matrix such that 
the I’s are brought together. Each block of l*s constitutes a 
component family and a machine cell. 

King [39] proposed a Rank order clustering (ROC) technique 
based on cluster analysis. This method interprets each entry in 
the matrix as a binary word and computes decimal equivalent of 
each row and column. First rows and then columns are sorted in 
decreasing order of the binary word values, until algorithm 
converges. 

Chandrasekharan and Rajagopalan [15] observed that ROC does - 
not provide block diagonal structure even in case of structured 
matrices where such a possibility exists. They found that ROC is 
dependent on the initial disposition of the incidence matrix and 
that it does not handle exceptional elements too well. They 


suggested MODROC. 



King and Nakornchai [40] observed that reading entries as 

binary words restricted the size of the problem to less than 47 

machines and 47 components because the largest integer 

48 

representation in many computers was (2 - 1). They proposed 

R0C2 which is computationally efficient. 

Chan and Milner [13] proposed direct clustering algorithm 
which is a poor version of ROC algorithm, except that it 
eliminates the sensitivity of the latter to the configuration of 
the initial matrix. This algorithm has same handicaps as the ROC 
algorithm and uses a limited binary comparison procedure for 
ranking. 

Graham [29] proposed linear weighting algorithm where the 
weights are increased linearly. The i th row is given a weighting 
of m - i + 1 where m is the total number of rows, and the 
priority ranking value is determined as the mean of the 
weightings of the non-zero entries. Ranking values calculated 
this way can be found and sorted very quickly and the requirement 
of a very large integer representation does not arise. The major 
disadvantages of. this linear weighting algorithm are the 
complicated and very confusing patterns of the intermediate 
results together with the difficulty in predicting the behavior 
of the procedure. 

McCormick [48] proposed Bond energy algorithm which reorders 
the rows and columns of the matrix for the purpose of moving 


numerically larger matrix elements together. A measure of 
closeness between two component vectors or machine vectors termed 
bond energy is defined. Method seeks to determine permutation of 
rows and columns in which the sum of products of adjacent 
elements is maximized. This is a restricted form of the quadratic 
assignment problem. 

Khator and Irani [38] proposed a heuristic procedure, the 
occupancy value method, for identifying clusters in a machine 
component matrix created from route card data. A unique feature 
of this method is that it progressively develops block 
diagonalization starting from the northwest corner of the matrix. 

Boe and Cheng [7] identified several deficiencies of 
clustering and array based sorting algorithms and suggested a 
clustering algorithm "close neighbour algorithm". In stage one 
the clustering of machines is performed and in stage two 
clustering of parts is done. The algorithm converges in one pass 
in most of the problems. 

Kusiak [41] proposed cluster identification algorithm (Cl) 
which identifies the cluster if it exists. He proposed three 
algorithms based on Cl to solve GT problem with bottleneck 


components and 

machines. 

The 

first 

algorithm 

solves an 

unconstrained GT 

problem. 

The 

second 

heuristic 

considers a 


constraint restricting the number of machines in each cell and 


identifies the bottlenecks. 



Tam [64] proposed a new similarity coefficient based on the 
operation sequence and used a clustering technique to obtain the 
component families. Such a coefficient, augmented with an 
advanced clustering technique, has been shown to improve 
productive effectiveness. 

2.6 Mathematical programming approaches 

Using mathematical approach of operations research (Linear 
programming, zero-one integer programming and dynamic 
programming), the GT cell production problem may be given concise 
mathematical formulation so that a solution can be found for a 
given constraint set. 

Gunasingh and Lashkari [30,31] proposed two zero-one integer 
programming formulations for the GT problem. The first 
formulation groups machines based on the compatibility of 
processing components, while the second formulation groups 
machines in order to minimize the cost of machine allocation and 
the cost of intercellular movement. 

Gunasingh and Lashkari [32] presented a non-linear 0-1 
integer programming formulation to simultaneously group machines 
and parts in cellular manufacturing systems based on tooling 
requirements of the parts, tools available on the machines and 
the processing times. The formulations take into account the 
limitations on the number of parts and machines in a group and 
the number of machine types available. 



Boctor [6] presented 0~1 integer programming based 
formulation for cell formation with the objective of minimizing 
the exceptional elements. This model allows the designer to 
control the cell sizes. 

Kusiak [43] developed a integer programming problem. He 
considered presence of alternative process plans and suggested a 
0-1 integer program based p-median model to determine the optimal 
process plan and assignments of parts to machine cells. 
Incorporating process plans will result in an improved quality of 
part families and the machine cells. 

Wei and Gaither [71] suggested an alternative formulation 
based on 0 - 1 integer programming for cell formation problem 
with the objective of minimizing the cost of manufacturing 
exceptional parts subject to constraint on the number of cells, 
cell size and the machine capacities. 

Kumar [66] have formulated a Joint grouping and loading 
problem as a multistage multi objective model. 

Shtub [58] showed that the generalized cell formation 
problem due to Kusiak [41] is equivalent to the generalized 
assignment problem. 

Choobineh [19] developed an integer programming problem to 
determine the assignments of component families to machine cells 
with the objective of minimizing the sum of the total production 
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cost and cost of purchasing new equipment. The model assumes that 
more than one component family be assigned to a cell. 

Srinivasan [60] proposed an assignment model for GT. The 
distance matrix for machines is solved as an assignment problem 
from which initial cells are identified. Component families are 
similarly identified and they are merged using a set of rules. 

2.7 Heuristic approaches 

Heuristic approaches have the ability to incorporate several 
practical considerations simultaneously. 

Waghodekar and Sahu [68] proposed an algorithm called 
machine component cell formation (MACE) to solve the grouping 
problem. The main feature of this heuristic is that it aims at 
minimizing the number of exceptional elements. 

Ballakur and Steudel [4] developed model which ' considers 
objectives such as minimizing exceptional elements, job-lateness 
and maximizing machine utilization imposing an upper limit on the 
number of machines in the cell. 

Askin and Subramaniam [2] developed a clustering algorithm 
that considers manufacturing costs like fixed and variable 
machining cost, setup cost, production cycle inventory cost, work 
in process inventory cost, material handling cost. The algorithm 
consists of three stages. In stage 1 parts are classified using a 
coding system. In stage 2 an attempt is made to develop a 
feasible grouping of parts based on the manufacturing costs. In 
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stage 3 the actual layout of machine cells is analyzed. 

Offodile [50] suggested a cell formation approach which aims 
at identifying parts with similar processing requirements. 

Harhalakis [35] described a twofold heuristic algorithm 
aimed at minimizing inter cell material movement. The first stage 
of the proposed heuristic is a bottom~up aggregation procedure to 
minimize the normalized inter cell traffic defined as the ratio 
of the inter cell movement between any two cells under 
consideration and the number of machines included in these cells. 
The second stage is aimed at improving the solution obtained from 
the first stage. 

Logendran [44] proposed a model which incorporates both 
intra cell and inter cell moves. Logendran [45] considered 
sequence of operations into account in the design of cellular 
manufacturing system. This model includes the impact of cells in 
evaluating inter cell moves. The model also aims at maintaining a 
targeted utilization of work station. 

Gupta and Seifoddini [33] proposed a methodology which 
incorporates production volume of various part types, routing 
sequence and unit operation times through an appropriately 
defined similarity coefficient. The first stage determines 
machine groups. The second stage assigns each part to one of the 


cells. 



Sule [63] outlined the importance of machine capacity 
planning in GT environment. Sule’s method determines the number 
of machines, their groupings and the amount of material transfer 
so that the cost of processing the components within the plant is 
minimized. 

2.8 Hybrid approaches 

The literature survey till now has shown that certain cell 
formation methods employs both heuristic and mathematical 
procedures in different stages. Such methods are called hybrid 
approaches. 

Steudel and Ballakur [4] proposed a two stage heuristic in 
which the first stage employs a dynamic programming procedure to 
generate an optimum chain of machines in which the sum of bonds 
between machines in the chain is maximized. The second stage is a 
heuristic, that partitions the chain to form machine cells 
subject to cell size restrictions. 

Nagi [49] addresses two problems (a) Process plan selection 
(formulated as a linear programming problem), (b) Cell formation, 
is solved through heuristic. The method also considers the 
sequence of operations and production data such as production 
volume, processing time, machine capacity. 

Choobineh [19] employed a two stage procedure. First employs 
a similarity measure based on part, operations and operation 
sequences and attempts to identify part families using a 



heuristic procedure. The second stage formulates an integer 
programming which specifies the type and number of machines in 
each cell and the assignment of part families to the cells. 

Askin and Chiu [1] proposed a mathematical model which addresses 
the problem of grouping of individual machines into cells and the 
routing of components to machines within cells. The mathematical 
programming formulation incorporates costs of inventory, machine 
depreciation, machine set-up and material handling costs. The 
formulation is then divided into two subproblems. The heuristic 
graph partitioning procedure has been proposed for each 
subproblem. The first subproblem assigns components to specific 
machines. The second sub problem groups machines into cells. They 
have also presented an approach for determining economic batch 
sizes in the group technology environment. 

2.9 Fuzzy theoretic approaches 

It was first advocated by Rajagopalan [54]. Chu and Hayya 
[20] observed that fuzzy approach reveals part families and also 
associates a degree of closeness between part and its part 
family. This information helps the designer in determining to 
which part family should be assigned as to balance workload among 
machine cells. 

2.10 Observations 

Following points were observed while doing literature 
survey. 

(1) Very less work has been done to integrate GT with materials 
requirements planning, just in time, computer integrated 
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manufacturing systems and with assembly operations. In my view 
such integration is essential to realize the benifits of group 
technology. 

(2) Very less efforts have been made to incorporate flexibility 
and flexibility measures in the design of group technology. 

(3) Group technology involves multiple objectives. Less efforts 
has been done to use mathematical techniques like goal 
programming in group technology. 

(4) Less work has been done in developing probabilistic 
approaches in the design of group technology. 

(5) There are many approaches available for design of group 
technology. The problem now arises is of finding which method is 
suitable for what type of problem. There is a need for comparison 
of group technology designs with this problem in mind. 

(6) Less work has been done to economically justify the methods 
of group technology. 

(7) Less work has been done to find out managerial implications 
on group technology. 

(8) Recent research is tending towards the application of 
artificial intelligence techniques and expert systems. 
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CHAPTER III 

SOME ALGORITHMS ON CELL FORMATION 

In this chapter some of the algorithms which are selected 
for subsequent comparison purpose are explained. These algorithms 
are both array based and non-array based heuristic algorithms. 
The same algorithms are coded for the purpose of comparison. The 
array based algorithms are 

(a) Rank order clustering (King [39]) 

(b) Linear weighting algorithm (Graham [29] ) 

(c) Direct clustering algorithm (Chan and Milner [13]) 

(d) Extended rank order clustering (King and Nakornchai [40]) 

(e) Bond energy algorithm (McCormick et. al. [48]) 
and the heuristic algorithms are 

(a) Cluster identification algorithm (Kusiak [41]) 

(b) Occupancy value method (Khator and Irani [38]) 

3.1 Rank order clustering (ROC) 

ROC is a well known clustering technique that attempts to 
create a block diagonal form by repeatedly reallocating the 
columns and rows of a machine part matrix according to binary 
values. King [39] shows that if the patterns of row entries are 
read as binary words then they can be ranked in reducing binary 
value order, using this principle ROC represents route card data 
as a binary matrix. Using a positional weighting technique for 
the *1' entries in the matrix, the rows and columns are 
alternately rearranged in order of decreasing rank. The result is 
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a diagonalization of the I’s into several clusters. If 
independent machine component groups do exist in the sample data 
provided, each machine will occur in only one cluster. Components 
will be uniquely assigned to any one of the clusters. Using this 
algorithm, the analyst can obtain a visual assessment of the 
machine groups and the associated families of parts 
simultaneously. With such an approach, a very valuable 
preliminary assessment of machines can be obtained because if a 
large number of machines are shared over several clusters, plans 
for cellular manufacture can be shelved at the out set. 


ALGORITHM 

STEP 1 : For each row of the machine~part incidence matrix, 


assign binary weight and calculate 

n 

For row i : a^^^ * 

k=l 

m 

For column j : Y * 

k=l 

m = no. of machines and n = no. of 

STEP 2 : Sort rows ' of the 

order of the corresponding decimal 

STEP 3 : Repeat the preceding two 


a decimal equivalent. 

n-k 

2 (3.1) 

m~k 

2 (3.2) 

parts. 

binary matrix in decreasing 
weights. 

steps for each column. 


STEP 4 : Repeat the preceding steps until the position of each 
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element in each row and column does not change. 

The steps of the algorithm are explained in Appendix through an 
illustration ( Appendix A.l ). 

3.2 Linear weighting algorithm 

Another approach for clustering data was proposed by Graham 
[29]. In this method the weights increases linearly according to 
the positions of the rows or columns. The i’th row is given a 
weighting of m - i + 1 whereas j’th column is given a weighting 
of n - J + 1. The priority ranking value is determined as the 
mean of the weightings of the non zero entries. Ranking values 
calculated this way can be found and sorted very quickly and the 
requirement of a very large integer representation does not 
arise. 


ALGORITHM 

STEP 1 : For each row of machine part incidence 

matrix assign linear weight. Calculate the weighted sum. 


For row i : * (n - k + 1) 


(3.3) 


k=l 


m 


For column j ; X! j * (m - k + 1 ) 
k=l 


(3.4) 


STEP 2 : Sort rows in decreasing order of corresponding 


weights. 
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STEP 3 : Repeat preceding steps until position of each element 

in each row and column does not change. 

The steps of algorithm are explained in Appendix through an 
illustration ( Appendix A. 2 ). 

3.3 Direct clustering algorithm 

The Direct clustering algorithm as proposed by chan and 
milner [13] is a poor version of the ROC algorithm, except that 
it eliminates the sensitivity of the latter to the configuration 
of the initial matrix. The number of positive cell entries K in 
each row and column is counted. The input machine component 
matrix is rearranged with columns in decreasing order of K but 
rows in increasing order of K. So, appropriate diagonalization is 
created and the machine component matrix is input in the same 
format always to the algorithm. 


ALGORITHM 

STEP 1 : Determine total number of Is in each row and column in 

machine part incidence matrix. 

n 

^ ^ (3.5) 

For row i : a^^ 

k = 1 


STEP 2 ; Sort each row in decreasing order corresponding to the 

total number of ones. 

STEP 3 : Sort each column in decreasing order corresponding to 
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the total number of ones. 

STEP 4 : Repeat the preceding steps until the position of each 

element and column does not change. 

The steps of algorithm are explained in Appendix through an 
illustration (Appendix A. 3) 


3.4 Bond energy analysis (BEA) 

McCormick et. al. [48] developed a matrix clustering 

technique which they call the bond energy algorithm. The matrix 

is applicable to any matrix in which non-negative integer values 

of an element in the matrix express a measure of the degree of 

association of the corresponding rows and column entities. The 

bond energy analysis algorithm attempts to identify and exhibit 

the interrelations within each cell and the associations among 

the clustered groups by means of total bond energy. A bond is 

claimed to exist between each pair of the neighboring rows and 

columns if they have positive cells in the machine-part matrix. 

BEA begins with an arbitrarily selected column (or rows). It then 

places that column with the greatest contribution to the total 

bond energy beside the assigned column (row). It repeats the same 

procedure for all the columns and rows. The BEA seeks to form a 

block diagonal form by maximizing the measure of effectiveness, 

which is defined as follows 
m n 


a. . . 

1-1. J 


.+ a 


i+1. j' 


(3.6) 
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ALGORITHM 

STEP 1 : Set j = 1. Select one of the columns arbitrarily. 

STEP 2 : Place each of the remaining n ~ j columns, one at a 
time, in each of j + 1 positions, and compute each columns 

contribution to the ME. Place the column that gives the largest 
incremental contribution to the ME in the best location. Increase 
j by 1 and repeat the preceding steps until j = n. 

STEP 3 : When all the columns have been placed, repeat the 

procedure for the rows. 

The steps of algorithm are explained in Appendix through an 
illustration ( Appendix A. 4 ). 

3.5 Extended Rank order clustering algorithm 

King and Nakornchai [40] observed that procedure of reading 

the entries as binary words presents some computational 

difficulties. Since the largest integer representation in most 
48 

computers is (2 - 1) or less, the maximum number of rows or 

columns that could be dealt with this way would be 47. They 
further observed that incidence matrices of the kind involved in 
group technology problems are usually very sparse, with densities 
unlikely to be higher than 5 to 10 %. This means that an 

elaborate system of linked list structures would in general be 
economical. 

King and Nakornchai [40] proposed new extended ROC algorithm 
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where the whole sorting procedure is reduced to that of shifting 
the orders of rows and columns. 

ALGORITHM 

STEP 1 : For each column, locate the rows (machines) with 

entries, move the rows with entries to the head of the row list, 
maintaining the previous order of the entries. 

STEP 2 : For each row, locate the columns with entries, move the 
columns with entries to the head of the column list, maintaining 
the previous order of the entries. 

STEP 3 : Repeat previous steps until there is no change in the 
matrix. 

The steps of algorithm are explained in Appendix through an 
illustration ( Appendix A. 5 ) 

2.6 Cluster identification algorithm 

In 1968 Iri [37] suggested masking technique. Starting from 
any row, mask all the columns which have an entry in this row, 
then proceed to make all rows which have entries in these 
columns. Repeat the process until the numbers of rows and columns 
stop increasing. These rows and columns constitute a block. 

In 1987 Kusiak and Chow [41] applied the concept presented 
in Iri [37] to develop the cluster identification algorithm 
(CIA). The cluster identification algorithm allows one to check 
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the existence of mutually separable clusters in a binary machine 
part incidence matrix provided they exist. 


ALGORITHM 

STEP 0 : Select iteration number k = 1. 

(k) 

STEP 1 : Select any row i of incidence matrix A and draw 

horizontal line h^ through it. 

STEP 2 : For each entry of 1 crossed by the horizontal line h^ 

draw a vertical line v.. 

J 

STEP 3 : For each entry of 1 crossed once by a vertical line 

Vj draw a horizontal line h^^. 

STEP 4 : Repeat steps 2 and 3 until there are no more crossed 

(k) 

once entries of 1 in A . All crossed twice entries of 1 in 
fk) 

A^ form machine cell MC-k and part family PF-k. 

(k) (k+1) 

STEP 5 : Transform the incidence matrix A into A by 

removing rows and columns corresponding to all the horizontal and 

vertical lines drawn in steps 1 through 4. 

STEP 6 : If matrix 0, where 0 denotes a matrix with all 

elements equal to zero, stop, otherwise set k = k + 1 and go to 

step 1. 

The steps of the algorithm are explained in Appendix through 
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an illustration ( Appendix A. 6 ). 

3.7 Occupancy value method 

Khator and Irani [38] proposed an heuristic approach to 
machine component grouping in cellular manufacturing 
applications. Using the assumption of independent machine 
component groups, a simple reordering of the listing of machines 
will decide the order in which the components get listed. The 
block diagonalization indicating clusters will appear. 

Occupancy value for a component can be calculated as: 

The route of a component, p^ consists of a set of machines 
I^. All components visiting one or more of these machines can be 
represented by a set J. Some of these components might be using 
additional machines, represented by a set I^. The machines in 
I^and ^2^^^ components in J represented a machine component sub 
matrix whose occupancy value is defined as 

OV = y M, / (m * n). (3.7) 

p jij J 

where m = number of rows of the sub matrix 

n = number of columns of the sub matrix. 

ALGORITHM 

STEP 1 : Develop an original machine component matrix 

consisting of ones and zeros. Form another initial machine 
component matrix identical to the original matrix. 
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STEP 2 : For each component j in the initial matrix find the 


total number of machines used 


STEP 3 : Scan the component list for selecting components with 
the minimum number of remaining machines (a) If only one 
component results from the above scan, go to step 5. (b) If there 
are ties in terms of more than one component having the same 
minimum number of operations go to step 4. 


STEP 4 : Calculate the Occupancy value for each of the 
components found in step 3. Select the component with the highest 
occupancy value. Ties for the highest occupancy value can be 
broken at random. 


STEP 5 : Enter the selected component and the machines used 

into the new matrix. 

STEP 6 : Update the initial machine component matrix for all 

the machines in the following manner. Cross out the rows 
containing machines selected in step 5. For each component 
remaining in the initial matrix decrease its Mj value by 1 if a 
crossed out machine was used by this component. 

STEP 7 : Enter any additional components in the new matrix if 

its M. value is zero. 

J 

STEP 8 : If all machines and components are entered in the new 

matrix, go to step 9, otherwise, go to step 3. 
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STEP 9 : Using cell entries from the original machine component 

matrix as input, compute the new matrix. 

The steps of the algorithm are explained in the appendix 
through an illustration (Appendix A. 7). 
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CHAPTER IV 

MEASURES OF COMPARISON 

Solution of practical problems without adequate 
generalizations and theoretical formulations is useful only in 
specific situations. Such results can not usually be applied to 
another problem. Group technology has suffered from this malady 
for a long time. When researchers and practitioners claim success 
with a certain method it is necessary to evaluate the result on 
an absolute quantitative scale. The practical use of such a 
measure is that subjective individual claims can be compared 
objectively. In case of block diagonal izat ion of zero-one 
matrices such a quantitative measure can be used as an objective 
function to be maximized. Such quantitative measures available in 
literature are : (a) grouping efficiency (GE) described by 

Chandrasekaran and Rajagopalan [14], (b) Measure of effectiveness 
as described by McCormick [48] and (c) Grouping efficiacy as 
described by Suresh kumar [18]. First quantitative measures, 
grouping efficiency, grouping efficiacy and Measure of 
effectiveness are described, followed by issues of bottleneck 
machines and exceptional parts. 

4.1 GROUPING EFFICIENCY 

The concept of grouping efficiency is developed to provide a 
quantitative standard on a rational scale for comparing different 
solutions to the same problem. 

Let {D >, where r = l,....k be the diagonal sub matrices 
r 
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obtained by associating the product group with the machine 
group (in the diagonalized matrix ). A perfect grouping 

implies that 

a = 1 if a € {D^} r = and a. = 0 if a, , s? {D }. 

iJ ij r ij ij r 

The diametrically opposite situation, therefore implies that 

a., = 0 if a., € {D} r = and a.. = 1 if a,. ^ {D }. 

iJ r ij ij r 

The efficiency scale is defined to cover these extreme cases, its 
value being set at 1 for the former case and zero for the latter. 

'Goodness* of grouping depends on two aspects: within group 
utilization and inter cell movement. From the matrix point of 
view, the concentration of non-zero elements in the diagonal sub 
matrices refers to utilization and the presence of such elements 
outside the diagonal sub matrices represents inter cell movement. 
A 'better grouping*, therefore, means an increase in the 
utilization, a decrease in the inter cell movement or both. 


Grouping efficiency can therefore be considered as a 
weighted average of two efficiencies and defined as 
follows. The number of non-zero elements in the diagonal blocks 
e^ is given by 


k r r 

EE ^ 

r=l 




(4.1) 


and the number of non-zero elements outside the diagonal blocks 
is Oq. Then 


= ®d ^ 


E 

r=l 


M 


N ) 
r 


( 


(4.2) 
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k 

The numerator of the expression for is the number of 
non-zero elements in the diagonal blocks and the denominator the 
total number of elements therein. Similarly, is the ratio of 
the number of zeros in the off-diagonal blocks to the total 
number of elements therein. 

The Grouping efficiency 7} can be expressed as the weighted 
average of and 7 )^. 

f] = + (1 - q)7]^ ( 4 . 4 ) 

where 0 ^ q 1 

7) satisfies the basic requirements of non dimensionality, 
non-negativity and zero to one range. The weighting factor q 
enables the analyst to alter the emphasis between utilization and 
inter cell movement, depending on the specific requirements of 
the given problem. 

4.2 LIMITATIONS OF GROUPING EFFICIENCY 

The weighted average of t }^ and gives the analyst enough 
freedom to decide the relative importance between inter cell 
movements and voids in diagonal blocks. Chandrasekharan and 
Rajagopalan contend that a value of q = 0.5 leads to the 
situation of equal weights to voids and exceptional elements. An 
analysis of the expression would reveal that this is not true. 
When large matrices are block diagonalized, for any number of 
blocks greater than two it can be observed that the first term in 
expression has a much smaller denominator compared to the second 
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term, whereas the numerators are more or less of the same order. 
As the matrix size increases, this disparity between the first 
and second terms widens and for large and/or sparse matrices the 
second term becomes less and less effective. In fact, if q and (1 
- q) are of the same order, the effect of inter cell movements 
(exceptional elements) is never reflected in the efficiency 
values in the case of large and sparse matrices. Therefore to 
have equal weights for voids and exceptional elements it is 
necessary to choose a very low value of q. Further, it is 
desirable to bring the two parts of the expression under a common 
denominator. 

4.2.1 EFFECT OF SPARSENESS OF THE MATRIX 

Choice of q would be more rational if it is linked with the 
size and sparsity of the matrix. As a result q is chosen as the 
density of an imaginary matrix that has an identical and perfect 
block diagonal structure as the given matrix. 

k 

q = y M N / mn. (4.5) 

r=l 

Placing value of q in the following equation 
n = q-ij^ + (1 - q)r)^. 
we get 

T) = 1 - (e- + e ) / mn. (4.6) 

0 V 

This simplified expression assumes that the perfect diagonal 
form has 100% efficiency which is reduced by the presence of 
voids and exceptional elements. 
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The following characteristics are important in the case of 
the modified expression for 7} after choosing q as discussed 
earlier. 

The main difference between this expression and the earlier 
expression for efficiency is that the present expression has the 
size of the matrix (mn) as its denominator, whereas in earlier 
expression 7}^and t \^ have two different denominators. It is also 
interesting to note that it is possible to assign relative 
weights for voids and exceptions in the present expression also, 
because it is in the present form that voids and exceptional 
elements have, an reality, equal weights. 

A major defect of grouping efficiency in both forms is its 

low discriminating power. It is obvious that (e^ is much 

less than the size of the matrix *mn*. Thus, for large matrices 

the term (e^ + e ) / mn will take values of the order of 10 which 
0 v 

leads to an efficiency close to one. 

Another serious defect lies in the definition of zero 
efficiency. Zero has been defined as the most degenerate case 
where e^ + e = mn. This implies that the diagonal blocks are 
full of zeros and off-diagonal blocks are full of ones. It can be 
seen that this definition is un-realistic from physical point of 
view. 

4.2.2 EFFECT OF SIZE OF THE MATRIX 

By accepting ’e* as denominator of the efficiency function it 
will have the form 

7) = 1 - (e^ + e^) / e. 
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“ (e - e - e ) / e. 

0 V 

= - e^) / e. (4.7) 

The main drawbacks of this function are as follows 
(1) Zero point is not properly defined, t] becomes 
zero when e , = e i.e. 


e 


d 


k 

% = 1/2 { I 
r = 


M N 
r r 




(4.8) 


This means that the zero efficiency is dictated by the fact 
that there are 50% voids in the diagonal blocks irrespective of 
the number of exceptional elements present. 

(2) This function takes negative values when e < e,. 

o V d 


4.3 GROUPING EFFICIACY 

For a large denominator the discriminating power of the 
efficiency function is poor and for a small denominator the 
function itself may become negative. 


The meaningful denominator between ’e* and *mn' is the 
’operational zone’ of the matrix. The operational zone of block 
diagonal form consists of the diagonal blocks usually consist of 
both ones and voids, and the exceptional elements consist of only 
ones. 

Operational zone = Voids + Ones in the diagonal blocks + Ones in 
of f -d i agona 1 blocks. 

= 0 + e , + e . 

V d 0 

= e + e. 

V 

Thus expression can be modified as 


(4.9) 



45 


T) = 1 - (e^ +6^) / (e^ + e) 

= (e^ + e - - e^) / (e^ + e) 

= (e - e„) / (e + e). (4.10) 

0 V 

This is grouping efficiacy (GEF). 

GEF = (1 - ^)/(l + $) (4.11) 

4^ = no. of exceptional elements 

Total no. of operations 
^ = no. of voids in diagonal blocks 

Total no. of operations 


4.4 MEASURE OF EFFECTIVENESS 

The measure of effectiveness (Measure of effectiveness) 
introduced in McCormick[48] is used to measure the quality of a 
solution. The measure is defined as follows: 


Measure of effectiveness 


n 


= I 


I I 

'1 1 


a. . 
1. J 


where a. . is the (i,J)th element 
1, J 

incidence matrix. 


* 


+ 

of 


Vi.j >■ 

the machine and part 


The Measure of effectiveness is specifically designed for 
data array clustering. This measure however is sensitive to the 
relative positioning of rows and columns of a matrix. The same 
grouping solutions with different arrangement of rows and columns 
may give different values of the measure. 
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4.5 EXCEPTIONAL PARTS AND BOTTLENECK MACHINES 

Most of the cell formation techniques employ part routeing 
information in the form of machine component 0-1 incidence 
matrix. A one as an entry in the matrix indicates that the 
corresponding machine is required to process corresponding part 
and a zero indicates otherwise. The matrix is rearranged to 
identify the machine cells and part families. 

Clustering of a binary incidence matrix may result in the 
following two categories of clusters: 

(1) Mutually separable clusters. 

(2) Partially separable clusters. 


0 1 3 2 4 5 

2 
4 

1 
3 

Figure 4.1 : Mutually Separable Cluster 
The above figure is Mutually separable cluster whereas the figure 
below is a partially separable cluster. 

0 
1 
2 

3 

4 

Figure 4.2 : Partially Separable Cluster 
There may be instances where a part may appear in more than 


1 2 3 4 5 
1 1 1 

1 1 

111 
1 1 
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one part family indicating that the processing of the part cannot 
be done in a single cell and thus an inter cell transfer would be 
required such parts are termed as exceptional parts. 

Matrix 2 cannot be separated into disjoint clusters because 
of part 5, which is to be machined in two cells, MC-1 and MC-2. 

Removing part 5 from matrix results in the decomposition of 

matrix into part families PF-1 = {1,2} and PF-2 = {3,4}. The two 

clusters are called the partially separable clusters. 

The presence of exceptional parts in the final solution of a 
cell formation problems may be due to the weakness of the 

algorithm or, in certain problem situations, it may be inevitable 
to avoid exceptional parts. 

Following actions to deal with exceptional parts can be 
taken. 

(1) It can be machined in one machine cell and transferred to 
the other machine cell by a material handling carrier. 

(2) It can be machined in a functional facility. 

(3) It can be sub-contracted. 

A bottleneck machine is a machine that does not allow 
decomposition of a machine part incidence matrix into sub 
matrices. 

Machine 3 in the figure 4.3 does not permit decomposition of 
that matrix into two machine cells and two part families. 
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Figure 4.3 : Effect of Bottleneck Machine 
A way to decompose matrix into two mutually separable sub 
matrices is to purchase an additional copy of machine 3. This 
will yield following figure. 
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Figure 4.4 : On Duplicating Bottleneck Machine 
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CHAPTER V 


IMPLEMENTATION, EXAMPLES, OBSERVATIONS AND EXPLANATIONS 


In order to compare the algorithms of cell formation 
discussed in chapter III several representative illustrations 
were choosen. Algorithms were coded using PASCAL and the system 
was implemented on HP 9000. Grouping efficiancy, grouping 
efficiacy, measure of effectiveness, number of bottleneck 
machines and number of exceptional parts were computed using the 
representative examples. Time for computation was recorded on HP 
system. Four representative examples which are selected for the 
discussions are given in section 4.1. Grouping efficiancy, 
grouping efficiacy, measure of effectiveness, improvement in 
measure of effectiveness, number of bottleneck machines and 
exceptional parts, and time are shown with each figure and at the 
same time are summarised in tables 1 to 4. 

5.1 REPRESENTATIVE EXAMPLES 



Figure 5.1 : Example 1 
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GE = 0.96 BM = 0 EP = 0 
Figure 5.2 : Output 
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ME =63 IME = 54 GEF 
by King’ s ROC Method 



GE = 0.96 BM = 0 EP = 0 ME = 63 IME = 54 GEF = 


0.92 


.92 


Figure 5.3 : Output by Graham’s LWA Method 
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GE = 0.66 BM = 0 EP = 10 ME = 28 IME = 19 GEF = 0.371 
Figure 5.4 : Output by DCA Method 



GE = 0.96 BM = 0 EP = 0 ME = 63 IME 


i . HAH'* 

» f. T.. 

A. « 


54 GEF = 0.92 


Figure 5.5 : Output by CIA Method 
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GE = 0.96 BM = 0 EP = 0 ME = 63 IME = 54 GEF = 0.92 
Figure 5.6 : Output by R0C2 Method 
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GE = 0.61 BM = 5 EP = 0 ME = 49 IME = 40 GEF = 0.348 
Figure 5.7 : Output by BEA Method 
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0.96 BM = 0 EP = 0 ME = 64 IME = 55 GEE 


Figure 5.8 : Output by OVM Method 
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Figure 5.9 : Example 2 
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GE = 1 BM = 0 EP = 0 ME 


68 IME = 46 GEF = 1 


Figure 5.10 


Output by ROC Method 



GE = 1 BM = 0 EP = 0 ME = 68 IME = 46 GEF = 1 
Figure 5.11 : Output by LWA Method 
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GE = 0.799 BM = 0 EP = 1 ME = 60 IME = 38 
Figure 5.12 : Output by DCA Method 
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= 1 BM = 0 EP = 0 ME = 68 IME = 46 GEF = 
Figure 5.13 ; Output by CIA Method 


GEF =0.6 
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GE = 1 BM = 0 EP = 0 ME = 68 IME = 46 GEF = 1 
Figure 5.14 : Output by R0C2 Method 
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GE = 0.81 BM = 0 EP = 0 ME = 62 IME = 40 GEF = 
Figure 5.15 : Output by BEA Method 
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1 BM = 0 EP = 0 ME = 68 IME = 46 GEE = 
Figure 5.16 : Output by OVM Method 
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Figure 5.17 : Example 3 
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GE = 0.75 BM = 1 EP = 1 ME = 58 IME = 26 GEF = 0.5 


Figure 5.18 : Output by ROC Method 
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GE = 0.63 BM = 1 EP = 3 ME = 42 IME = 10 GEF = 0.3 
Figure 5.19 : Output by LWA Method 
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GE = 0.6 BM = 1 EP = 14 ME = 28 IME = 0 GEF = 0.24 
Figure 5.20 : Output by DCA Method 
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GE = 0.67 BM = 0 EP = 0 ME = 50 IME = 18 GEF = 0.33 
Figure 5.21 : Output by CIA Method 
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Output by R0C2 Method 
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GE = 0.74 BM = 0 EP = 1 ME = 57 IME = 25 GEF - 0.48 
Figure 5.23 : Output by OVM Method 
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Figure 5.24 : Example 4 
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GE = 0.65 BM = 1 EP = 1 ME = 113 IME = 57 GEF = 0.37 
Figure 5.25 : Output by ROC Method 
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GE = 0.58 BM = 1 EP = 14 ME = 80 IME = 24 GEF = 0.29 


Figure 5.26 : Output by LWA Method 
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GE = 0.56 BM = 0 EP = 22 ME = 59 IME = 3 GEF - 0.262 
Figure 5.27 : Output by DCA Method 
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OUTPUT BY KING’S ROC2 METHOD 
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GE = 0.68 BM = 2 EP = 6 ME = 113 IME = 57 GEF = 0.40 


Figure 5.28 : Output by R0C2 Method 
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GE = 0.57 BM = 7 EP = 14 IME - 55 GEF 0.276 
5.29 : Output by BEA Method 


Figure 
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GE = 0.70 BM = 1 EP = 4 ME = 98 IME = 42 GEF = 0.42 


Figure 5.30 ; Output by OVM Method 
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TABLE 1 ; COMPARISON OF VARIOUS ALGORITHMS 




1 BM = Bottleneck machines 

2 EP = Exceptional parts 

3 GE = Grouping efficiancy 

4 ME = Measure of effectiveness 

5 IME = Improvement in measure of effectiveness 

6 GEE = Grouping efficiacy 
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TABLE 3 : COMPARISON OF VARIOUS ALGORITHMS 


PROBLEM 

SIZE 


ROC LWA DCA CIA ROC2 BEA OVM 
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TABLE 4 : COMPARISON OF VARIOUS ALGORITHMS 
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4.2 Observations 


The following observations were made from the results 

(1) Rank order clustering [39], extended rank order clustering 
[40] and occupancy value method [38] gave nearly equal grouping 
efficiencies. 

(2) Cluster identification method [41] was successful in giving 
good grouping efficiency for some example, comparable with kings 
[39] and Occupancy value method, but it failed altogether for 
some examples. Cluster identification algorithm failed to 
identify clusters when machine part matrix was dense, and had 
more number of bottleneck machines and exceptional parts. 

(3) The grouping efficiency of linear weighting algorithm [29] 
followed that of Rank order clustering and Occupancy value 
method, which was followed by Direct clustering algorithm [13] 
and then by bond energy algorithm [48], 

(4) When a large matrix of 16 * 43 machines and parts was fed to 
the program. Rank order clustering [39] failed to group, because 
of integer overflow, whereas other methods gave answers. Thus 
Rank order clustering is not useful for larger matrices. 

(5) For some examples Rank order clustering converged to solution 
which were not optimal but such occurrence were very few as 
compared to Linear weighting algorithm. Direct clustering 
algorithm and Bond energy algorithm. 
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(6) In Linear weighting algorithm the ranking values were sorted 

quickly and requirement of a very large Integer representation 

does not arise. 

(7) Linear weighting algorithm produced confusing patterns of the 
intermediate results together with the difficulty in predicting 
the behavior of the procedure. 

(8) Linear weighting algorithm and Direct clustering algorithm 
were not able to produce exact diagonal matrices in most of the 
cases. The examination of resulting matrices may sometimes 
suggest one or two column or row interchanges to achieve a 
desired result. 

(9) Bond energy algorithm gave solution matrix with a 
chequerboard structure with overlapping blocks. 

(10) Bond energy algorithm gave more number of bottleneck 
machines than any other method. 

(11) Direct clustering algorithm gave maximum number of 
exceptional parts followed by Bond energy algorithm. 

(12) Rank order clustering, Extended Rank order clustering and 
Occupancy value method gave highest measure of effectiveness, it 
was followed by Cluster identification algorithm. Linear 
weighting algorithm, Direct clustering algorithm and Bond energy 
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algorithm gave relatively poorer measure of effectiveness. 

(13) Cluster identification algorithm cannot deal with problems 
containing bottleneck parts. 


(14) When we compare the CPU time Direct clustering algorithm, 
Bond energy algorithm, Linear weighting algorithm and Cluster 
identification algorithm takes less time as they involve less or 
no computations whereas Rank order clustering. Extended rank 
order clustering and Occupancy value method takes longer time as 
computations are involved. 

4.3 Explanation 

Rank order clustering [39] uses exponential weights and 
binary ranking. The reason why Rank order clustering gives good 
solution is that it assigns positional weights and that weight is 
of substantial order. There are less chances of ties in ranking. 
Mostly the tie will be only in case of exactly similar rows or 
columns. Thus the possibility of alternate optimum solution 
possibility is less. The positional weight gives justice to all 
the elements and binary ranking becomes discriminatory. 

In some cases Rank order clustering does not give optimal 
results. This may be because of ties in columns and rows which 
allows alternate and better solution if rows and columns are 
rearranged. Rank order clustering method may not produce 
consistent results, when rows or columns in the initial matrix 
interchanged. Thus Rank order clustering depends on initial 


are 



74 


algorithm gave relatively poorer measure of effectiveness. 

(13) Cluster identification algorithm cannot deal with problems 
containing bottleneck parts. 

(14) When we compare the CPU time Direct clustering algorithm. 
Bond energy algorithm, Linear weighting algorithm and Cluster 
identification algorithm takes less time as they involve less or 
no computations whereas Rank order clustering. Extended rank 
order clustering and Occupancy value method takes longer time as 
computations are involved. 

4.3 Explanation 

Rank order clustering [39] uses exponential weights and 
binary ranking. The reason why Rank order clustering gives good 
solution is that it assigns positional weights and that weight is 
of substantial order. There are less chances of ties in ranking. 
Mostly the tie will be only in case of exactly similar rows or 
columns. Thus the possibility of alternate optimum solution 
possibility is less. The positional weight gives justice to all 
the elements and binary ranking becomes discriminatory. 

In some cases Rank order clustering does not give optimal 
results. This may be because of ties in columns and rows which 
allows alternate and better solution if rows and columns are 
rearranged. Rank order clustering method may not produce 
consistent results, when rows or columns in the initial matrix 
interchanged. Thus Rank order clustering depends on initial 


are 
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disposition of the patrls. This occurs because, by changing rous 
and columns we are changing the positional weights which leads to 

different answers. 


The occurrence of exceptional elements is expected but Rank 
order clustering solution is disrupted, due to the method adopted 
for ranking. It relies on pair wise comparisons of cell entries 
in the leftmost column (when ranking rows) and topmost row (when 
ranking columns). So if the positional occurrence of these 
elements is such that they influence the ranking, poor cluster 
formation will result. Bottleneck machines are used by a large 
number of components. Since these components can be expected to 
be dispersed over more than one cluster, such machine must appear 
in more than one row in the matrix. Otherwise the ranking 
procedure creates large, dispersed clusters with many machines 
and components contained in them. Rank order clustering algorithm 
will work more effectively after bottleneck machines and 
exceptional parts are identified and suppressed after visual 
analysis of the initial matrix solution. Such prior assumptions 
bias the solution, especially as the algorithm must indicate 
exceptions and bottleneck machines, not rely on their temporary 
suppression to be effective. The reason why sometimes Rank order 
clustering fails is because king [39] assumes that bottleneck 
machines can be freely duplicated if they are required in several 
clusters. He ignores the fact that some clusters would need to be 
merged to optimize the utilization of these machines. 


Rank order clustering fails for larger 


The reason why 
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matrices is because of the procedure of reading the entries as 
binary words which presents some computational difficulties, 
since the largest integer representation in most computers is 
(2 - 1 ) or less, the maximum number of rows or columns that 
could be dealt with this way would be 47. Also the storage of the 
Incidence matrix as a two dimensional array puts a severe limit 
on the problem that can be tackled. A moderate problem with 50 
machines and 2000 components, together with the program would 
require core storage in excess of 120 K. Secondly because the 
sorting procedure has a complexity of cubic order, efficient 
implementation is not possible for very large problems. Extended 
Rank order clustering [40] was proposed to overcome this 
limitation. It uses the observation that GT problems are usually 
very sparse, with densities unlikely to be higher than 5 to 10 % 

. Extended Rank order clustering uses elaborate system of linked 
list structures and is very successful. Thus extended Rank order 
clustering gives answers when Rank order clustering fails. 
Further Rank order clustering and extended Rank order clustering 
uses array manipulations which reduces their efficiencies. 
However Rank order clustering gave good results comparable with 
Occupancy value method [38]. 

Occupancy value method gives good results because 

(a) It creates compact clusters without any preliminary 
assumptions or visual identification of exceptional elements or 
bottleneck machines. 

(b) It builds up clusters along the diagonal by using small 
of the larger original matrix. This allows the 


selected sections 
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analyst flexibility to re-iterate when ever dispersion of cell 
entries outside a cluster is observed, due to the occurrence of 
bottleneck machines. So, given any seed or starting machine, its 
components will indicate the other machines that will occur in 
the cluster, only a good choice of seed machine is required. 

(c) It uses the machine component matrix only as a means for data 
representation and not for array manipulations like the Rank 
order clustering, Linear weighting algorithm, and Direct 
clustering algorithm. 

(d) It uses the occupancy value to delay the entry of bottleneck 
machines and components with a large number of machines in their 
routes into the matrix. Otherwise this will create large 
clusters, useless for transforming into cells. 

(e) It simplifies the identification of both exceptional elements 
and bottleneck machines by grouping correctly all the machines 
which would occur in only one cluster. 

In Cl (Kusiak [41]) the results obtained are good, but for 
many matrices with more density and more exceptional parts and 
bottleneck machines, the algorithm fails. This is because of the 
nature of algorithm which uses line drawing procedures. If matrix 
is dense all the rows and columns gets crossed resulting in one 
cell, the same as original matrix. 

In Linear weighting algorithm (Graham [29]) the weighting of 
positions is linear. The weights increases linearly. This method 
can solve large matrices because requirement of a very large 
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obtained are Inferior to Rank order clustering and Occupancy 
value method because it uses small positional weight for rows and 
columns. Hence ranking may involve ties. This leads to confusing 
patterns of the intermediate results together with the difficulty 
in predicting the behavior of the procedure. This method produces 
results inferior to Rank order clustering and Occupancy value 
method but better than Direct clustering algorithm and Bond 
energy algorithm. 

In Direct clustering algorithm we do not use positional 
weighting procedure. We sum the entries of ones for each row and 
column and sort it according to their sum value. Direct 
clustering algorithm does not work properly because of ties in 
ranking which is very predominant. Adjustment of rows and columns 
may yield better results. It has all the difficulties of Rank 
order clustering except the storage. 

The number of entries in the rows and columns which will 
increase in the matrices analyzed later, creating cluster 
dispersion similar to the influence of the bottleneck machines 
and exceptional elements. This method allows more flexibility in 
the size of problem. Furthermore, the sensitivity of the Rank 
order clustering algorithm to the initial matrix is eradicated 
because Direct clustering algorithm initiates the procedure by 
counting the number of positive cells instead of depending on 

intuition. 

Bond energy algorlth. (48) Is applicable to problems of any 



79 


size because Bond energy algorithm has nothing to do with 
calculating the binary values. However the first step (select row 
or column) is determined by intuition. Many possible solutions 
can be generated, i.e the solution depends on the initial row or 
column selected for starting the process. Bond energy algorithm 
gives chequer board structure because the main aim of the 
algorithm is to improve measure of effectiveness and not to 
obtain a diagonal structure. The algorithms groups ones but they 
are not grouped in diagonal manner only the close proximity of 
the ones is enough to converge. 
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CHAPTER VI . 


CONCLUSIONS 


6.1 Conclusion 

In this thesis role of group technology in booasting 
the productivity and its implications on the dynamic business 
environment is discussed. The concept of cell formation and 
division of system into subsystem has been discussed at 
appropriate places. A comprehensive literature survey has been 
presented and observations has been made on literature available. 
It has been stressed that although vast literature is available 
on group technology, comparison of different techniques on common 
scales has not been made. Many algorithms are explained in 
details giving the examples. These algorithms are both array type 
and heuristic algorithms. The algorithms are coded and results 
for cell formation are obtained therefrom. Four representative 
examples are choosen and algorithms are compared on basis of 
various measures like grouping efficiency, grouping efficiacy, 
measure of effectiveness, improvement in measure of 

effectiveness, time, number of bottleneck machines and 

exceptional parts. This measures of comparison are discussed in a 
separate chapter. The results, tables and graphs are used as an 
aid for comparison. It has been found that rank order clustering, 
extended rank order clustering, occupancy value method and 
cluster identification method gives similar results which are 
better than direct clustering algorithm, linear weighting 
and' bond energy algorithm. The observations. 


algorithm. 
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explanations and why different algorithms behave in a particular 
way, are given. Thesis ends with an observation that much can be 
done in this field by comparing more and more algorithms to find 
the best and suitable technique. 

5.2 Scope for further research 

Despite numerous economic benifits and operational 
advantages offered by the group technology, its real potential 
has not been fully explored. 

Alternative approaches using recent techniques such as 
nueral nets, genetic algorithms, tabu search etc., in obtaining 
machine cells can be examined. Work can be done to integrate cell 
formation techniques with other functions of manufacturing 
organization such as design, process planning, production 
control, materials requirements planning, computer integrated 
manufacturing systems etc. 

Various types of flexibilities can be incorporated in 
design of cellular manufacturing system. Use of goal programming 
for solving multiple objectives can be analysed. Probabilistic 
approaches, fuzzy theoretic approaches, polyhedral dynamic 
approaches, and expert system can be used in group technology. 
There is a need to exploit the power of simulation by simulating 
potential cellular manufacturing designs, before a decision is 
made on the selection of a particular design. Economic 
justification can be incorporated into group technology and more 
managerial issues can be analysed. 
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APPENDIX A 

A.l 

This section gives the illustration for the Kings Rank order 
clustering [39] method. The original matrix given is as follows 
(Figure 1) on which the computations are done to obtain the final 
matrix. Further computions are done on the following figures. 

Binary 6 5 4 3 2 1 0 decimal binary 

weight 2 2 2 2 2 2 2 weight ranking 

0 ^ 2 3 4 ^ ^ 7 

1 0 1 0 1 1 1 0 46 4 

21010000 80 2 
3 1 0 1 0 0 0 1 81 1 

4 0 1 0 1 0 1 0 42 5 

5 [l 0 0 0 0 0 1 I 65 3 

Figure 1 : Initial Matrix and ROC Calculations 
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Figure 2 : Calculation for ROC Method 
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Figure 3 ; Final Matrix by ROC Method 

Figure 3 is the final matrix obtained by the method. 
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A . 2 

This section illustrates the grahams linear weighting 
algorithm, the same illustration is used for solving by this 
method. 
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Figure 4 : Calculation of LWA Method 
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Figure 5 : Calculation of LWA Method 
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Figur'c 6 t Calculation of LV/A Method 
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Figure 8 : Calculation of LWA Method 
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Figure 9 : Block Diagonal Structure by LWA Method 


This section illustrates chan and Milners. Direct Clustering 
Algorithm [13] with the same example. 
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Figure 10 : Calculation by DCA Method 
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Fi^niic 11 : Calculation by DCA Method 

^ '/ 3 4 6 2 5 

2 1 0 1 0 0 0 0 

5 1 1* 0 0 0 0 0 

3 1 1 1 0 0 0 0 

4 0 0 0 1 1 1 0 

1 0 0 0 1 1 1 1 

Igur.r 12 : Block Diagonal Structure by DCA Method 


JTH. w 

Thif; section lllustratres the Cluster Identification 
m,™ o< Kuslak (411. Selecting the first rou as the selected 
row in the figure 1 «e apply line drawing procedure to get 

in-alal sot as .achlnes (1. 4) and the part set as the 

,,a,ts (;-.4,S,C>. further choosing the .achlne 2 and applying the 
line drawing procedure we get .achlne set (2,3.5) and the part 
set as (1.3.5). the final .atrlx obtained Is 
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Figure 13 : Block Diagonal Structure by CIA Method 

A. 5 

This lUH'tlon illustrates the Kings R0C2 algorithm [40]. The 
algoi Hhm Is applied on figure 1. 
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14 : Calculation of R0C2 Method 
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For row no. 


column list 

5 1 2 3 4 5 6 7 

42461357 
32465137 
21724653 
1 1 3 7 2 4 6 5 


New column order 1 372465 

Figure 15 : Calculation by R0C2 Method 


Th<‘ f inal matrix obtained is as follows 
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Figure 16 : Block Diagonal Structure by R0C2 Method 


1 3 __ 7 2 4 6 5 
1 1 1 0 0 0 0 
1 1 0 0 0 0 0 
1 0 1 0 0 0 0 
0 0 0 1 1 1 1 
0 0 0 1 1 1 0 


A. 6 

Thin section illustrates the Bond energy method of 
McC'oi mi ek 1 48 ) . The computations are done on figure 1. 

Step 1 : Set J = 1. Select column 2. 

Step 7. : Place each of the remaining columns in each of the j +1 
positions. The contributions to the ME value of column 2 is 


computed next. 
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FiHun? 17 : Calculation by BEA Method 


Column 4 i*; [ilatod in the J + 1 th position. 

Ing thf same pi ocodure for the remaining columns leads to 
the following column order {2,4,16,6,1,3,7}. 

Step 3 : Repeat lag the preceding steps for rows results on the 
row ord»“i {1,3, 2, 4, 5}. The final matrix obtained is 


Flguie IH 


C) 
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?. J 6 5 13 7 

11110 0 0 
0 0 0 0 1 1 1 

0 0 0 0 1 1 0 

1 1 1 0 0 0 0 

0 0 0 0 1 0 1 

mock Diagonal Structure by BEA Method 


A. 7 

This section illustrates the method of solving a problem with 
the heflp of occupancy value method. The original matrix is given 

in i igure 1. 
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rn;HATi(iN i 

"■'•■P ■> ■ iui Identical initial matrix [I] for use in all 

*;ul>::*‘<jiifnt ntejt;;. 


Stop ?. 

: i"ca 

each component 

J in [I] 

scan 

t ot a I ruimb<‘i 
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Stfp : Scan Mj values in component list for selecting 
cnmptincnt wltli tiu? minimum number of machines. Only one is found 
tliat it; ( ompuiifiit 5. 

Stfj> 5 ; Comiinncnt S nad machine 1 enters the new matrix. 

Stop 6 : Ttio initial matrix is now updated. Set all a values in 

1 j 

the rows; coi r etipomding to the entering machines zero and compute 
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Figure 19 ; Calculation by OVM Method 
Iteration 2 

Step 3 : Minimum Mj value is 1. 

step 4 : Since ell three exponents require the sane nachlne 4 
occupancy value will be equal for all. 


Step 5 : Enter any one 


of the component (say 2) and machine 4 
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into new matrix. 

Step 6 : Update the matrix and find 
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Figure 20 : Calculation by OVM Method 


ITERATION 3 

Step 3 : Minimum Mj found is 2. There are ties so we go to step 

4. 


Step 4 : We find occupancy value for component 3 


Ij = {2,3} I 2 = {5> J = {1.3,7} 

0 
2 
3 
5 

M. 3 2 2 

J 

OV = S My (m* n) 

= 7/9 

= 0.777 


1 3 7 
110 
1 1 1 
1 0 1 


Occupancy value for the component 7 will be same. 
We enter component 3 and machines 2 and 3. 





98 


0 

1 

2 

3 

4 

5 

6 

7 

1 

0 

0 

0 

0 

0 

0 

0 

2 

0 

0 

0 

0 

0 

0 

0 

3 

0 

0 

0 

0 

0 

0 

0 

4 

0 

0 

0 

0 

0 

0 

0 

5 

1 

0 

0 

0 

0 

0 

1 


Figure 21 : Calculation by OVM Method 
Next we enter component 1 and 7 and machine 5, which yields 
full zero matrix. Final new matrix obtained is 
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Figure 22 : Block Diagonal Structure by OVM Method 


A. 8 

This section illustrates the method of calculating the 

grouping efficiency and grouping efficiacy. 

Grouping efficiency ; 

7) = qT7j + (1 - q)V2 

V = e y I M N = 14 / 17 = 0.8235 

*1 d ^ r r 

7)2 = 1 - Sq / (mn - Z Mj. N^) = 1 
Taking q = 0.5 we get t) = 0.911. 
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Figure 21 : Calculation by OVM Method 
Next we enter component 1 and 7 and machine 5, which yields 
full zero matrix. Final new matrix obtained is 
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Figure 22 : Block Diagonal Structure by OVM Method 


A. 8 

This section illustrates the method of calculating the 

grouping efficiency and grouping efficiacy. 

Grouping efficiency : 

7} = + (1 - q)'i?2 

7) = ey E M N = 14 / 17 = 0.8235 

1 d r r 

7)2 = 1 - Sq / (mn - E Mr ^r^ ^ ^ 

Taking q = 0.5 we get t) = 0.911. 
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Grouping efficiacy 




4 > 


GEF = (1 - \&) / (1 + ,^) 


No. of exceptional elements 
Totalno. of operations 


0 

14 


0 . 


No. of voids in diagonal blocks 
Totalno. of operations 



0.2143. 


GEF = 1 / 1.2143 = 0.8235. 



